High energy electroweak DVCS processes in neutrino-nucleon scattering 
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In this work we estimate the differential and total cross sections for the high energy deeply virtual 
Compton scattering in the weak sector. In the weak neutral sector one considers neutrino scattering 
off an unpolarized proton target through the exchange of Z°. We numerically compute the process 
Z*p — * yp within the QCD color dipole formalism, which successfully describes the current high 
energy electromagnetic DVCS experimental data. We also discuss possible applications for the weak 
charged sector and perform predictions for scattering on nuclear targets. 
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I. INTRODUCTION 

The cross section of hard scattering processes can be 
written as the convolution of parton distributions (PDFs) 
and the cross sections of hard subprocesses computed at 
parton level using perturbative QCD. The usual PDFs 
are the diagonal element of an operator in the Wilson op- 
erator product expansion, which can be extracted from 
data for inclusive processes. On the other hand, there are 
a set of exclusive reactions which are described by the 
off-diagonal elements of the density matrix, where the 
momentum, helicity or charge of the outgoing target are 
not the same as those of the correspondent ingoing par- 
ticle. Examples of such reactions are the virtual photon 
Compton scattering and the vector meson electroproduc- 
tion. In these cases, the difference with the inclusive case 
is the longitudinal components of the incoming and out- 
going proton momentum. Precision data are becoming 
available for hard scattering processes whose description 
requires knowledge of these generalized parton distribu- 
tions (GPDs). 

In recent years, significant effort was made to access 
GPDs through the measurement of hard exclusive elec- 
troproduction processes. The simplest process in this 
respect is the electromagnetic deeply virtual Compton 
scattering (DVCS), that is an electron scatters off a nu- 
cleon producing a real photon in the final state j*p — > jp. 
On the other hand, studies on exclusive processes medi- 
ated by weak interactions have increased motivated by 
experimental proposals for neutrino factories and high in- 
tensity neutrino beams facilities. Using a neutrino beam 
instead of an electromagnetic one the DVCS process into 
the weak sector can be investigated. For the neutral weak 
current one has the process Z°p —> jp, whereas for the 
charged current case occurs the process W^p — > jN. In 
these cases, one expects sensitivity to a different flavor 
decomposition of GPDs and also due to the presence of 
the axial part of the V-A interaction they are sensitive 
to a different set of GPDs. Therefore, neutrino scat- 
tering on nucleons in the deeply virtual kinematics can 
shed light on the hadronic interactions mediated by elec- 
troweak currents. 



process in case of weak interactions. In Ref. [l[ the neu- 
trino scattering off protons in the deeply virtual kinemat- 
ics, which describes under a unified formalism elastic and 
deep inelastic neutrino scattering, was first studied via an 
extension of the notion of DVCS to the case of a neutral 
current exchange. Afterwards, in Ref. Q the investi- 
gation was extended for both charged and the neutral 
electroweak sectors, where the structure of the leading 
twist amplitudes of the process has been discussed. In a 
very updated investigation the DVCS in weak sector 
is considered in the generalized Bjorkcn limit and calcu- 
lated at the leading twist within the framework of the 
nonlocal light-cone expansion via coordinate space QCD 
string operators. The cross sections for neutrino scat- 
tering off the nucleon, relevant for future high-intensity 
neutrino beam facilities are estimated. 

Besides the GPDs formalism for exclusive process, at 
high energies a simple and intuitive approach has at- 
tracted a wide interest. It is the color dipole formal- 
ism [ij, where the interaction is described in the target 
rest frame and relies on the Fock states expansion of the 
projectile and their further scattering with the target. 
Specifically, the exclusive process (7, Z°, W )p — > Ep 
is described by the interaction of qq pairs (color dipoles), 
in which the virtual boson fluctuate into, with the nu- 
cleon. The scattering amplitude is given by the convolu- 
tion of the dipole-nucleon cross section with the overlap 
of the virtual incoming boson and outgoing exclusive (E) 
final state wave functions. This approach provides a very 
good description of the data on yp inclusive production, 
77 processes, diffractive deep inelastic and vector meson 
production (For a recent review, see [f|). In particular, 
the electromagnetic DVCS cross section is nicely repro- 
duced in several implementations of the dipole cross sec- 
tion at low cc @, 0] ■ In particular, the ability to predict 
cross sections at very high energies depends critically on 
the reliability of extrapolations from current measure- 
ments into regions of lower virtualities of the elementary 
scattering process, and much smaller values of the frac- 
tional momentum x carried by the parton constituents of 
the hadrons. The color dipole approach is quite robust 
for doing this job. 



Recently, some sound works have addressed the DVCS 



In this letter we study the application of color dipole 
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formalism to the DVCS process in electroweak sector. 
For real photon in the final state, the QCD calculation 
becomes not trivial and saturation models can provide 
a suitable scheme for dealing with low momenta scales 
appearing in the problem. The present calculation is 
strongly motivated by recent studies in Refs. [H, 0, 0|, 
where the leading twist amplitudes for deeply virtual 
neutrino scattering (DVNS) have been discussed. Those 
investigations have considered the generalized Bjorken 
limit within the framework of the nonlocal light-cone ex- 
pansion, which involves GPDs only for quarks. On the 
other hand, at higher energies the processes induced by 
gluons are increasingly important. Therefore, the stud- 
ies presented here are complementary to those in Refs. 
P, [3, El- In next section, we make a short summary 
on the color dipole approach applied to electromagnetic 
DVCS and show its success in describing the high energy 
accelerator data. In Sec. 3, we extend the formalism for 
the electroweak sector, considering for simplicity the neu- 
tral current case. Prediction are done for the differential 
cross section da/dt and total cross sections considering a 
proton target. In addition, phenomenology considering 
nuclear targets is discussed. Finally, in last section we 
summarize the results. 



II. HIGH ENERGY ELECTROMAGNETIC 
DVCS AND THE COLOR DIPOLE APPROACH 

Let us introduce the main formulas concerning the 
color dipole picture applied to DVCS processes. In the 
dipole model deep inelastic scattering is viewed as 
the interaction of a colour dipole, i.e. mostly a quark- 
antiquark pair, with the proton. The transverse size of 
the pair is denoted by r and a quark carries a fraction z 
of the photon's light-cone momentum. In the proton rest 
frame, the dipole lifetime is much longer than the life- 
time of its interaction with the target proton. Therefore, 
the elastic 7*2? scattering is assumed to proceed in three 
stages: first the incoming virtual photon fluctuates into 
a quark-antiquark pair, then the qq pair scatters elasti- 
cally on the proton, and finally the qq pair recombincs 
to form a virtual photon. The amplitude for the elastic 
process 7*p — > j*p, A 1 p (x, Q, A), is simply the product 
of amplitudes of these three subprocesses integrated over 
the dipole variables r and z: 

A 7 " p (x,Q,A) = rV fd 2 r [ dz% n (r, Z> Q) 
f h,h J J ° 
x A qq {x,r,A)<$ h - h (r,z,Q), (1) 

where ^^(r, z, Q) denotes the amplitude for the incom- 
ing virtual photon to fluctuate into a quark-antiquark 
dipole with hclicitics h and h and flavour /. The quantity 
A q q{x, r, A) is the elementary amplitude for the scatter- 
ing of a dipole of size r on the proton, A denotes the 
transverse momentum lost by the outgoing proton (with 
t = —A 2 ), x is the Bjorken variable and Q 2 is the photon 



virtuality. The elementary elastic amplitude A qq can be 
related to the S'-matrix element S(x, r, b) for the scatter- 
ing of a dipole of size r at impact parameter b [H : 



A qq (x,r,A) 



Jd 2 be- ib *A qq (x,r, b) 



d 2 b 



,-ib-A 



2[1-S(x )r> b)].(2) 



This corresponds to the notion of impact parameter when 
the dipole size is small compared to the size of the proton. 
The integration over b of the S-matrix element motivates 
the definition of the qq-p differential cross section as 



da„ 



d 2 b 



2[1- Re S(x,r, b)}. 



(3) 



The amplitude for production of the exclusive final 
state such as a real photon in DVCS, is given by 

A^ p ^ p {x,Q,A) = [d 2 r f dz(%%,) T A qq (x,r,A) 



= i d 2 r 



( dz J d 2 b (^;* 



7-Jt e 



-ib-A d<J qq 



d 2 b 



(4) 



where (^*'^>)t denotes the overlap of the vitual incoming 
photon and outgoing real photon wave functions. For 
DVCS, the amplitude involves a sum over quark flavours. 
Here, one assumes that the size of the quark-antiquark 
pair is much smaller than the size of the proton. After 
considerations about the non-forward wave functions we 
finally obtain 0], 



A 



7 p^7p 



(x,Q,A) = i /d a r / dz jd 2 b{%y r ) T 



-i[b-(l- Z )r]-A dg gg 
X 6 d 2 6' 

The elastic diffractive cross section is then given by 



dlj 7 p^jp 
dt 



1 

16^ 



A^ p ^ p (x,Q,A) 



(5) 



(6) 



For one has real photons in final state, only the trans- 
versely polarized overlap function contributes to the cross 
section. Summed over the quark helicities, for a given 
quark flavour / it is given by 



(*;*7*)^ 



2ir 2 



{[z 2 + z 2 ] eK^erjmfKximfr) 



(7) 



+ m 2 K (er)Ko(m f r)} , 

where we have defined the quantities e 2 = zzQ 2 
and z = (1 — z). 

Before discussing the specific model for the dipole-cross 
section some corrections to the exclusive processes should 
be shortly addressed. For DVCS one should use the off- 
diagonal (or generalized) gluon distribution, since the ex- 
changed gluons carry different fractions x and x' of the 
proton's (light-cone) momentum. The skewed effect can 
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FIG. 1: (Color online) Differential cross section da/dt for 
DVCS on a proton as a function \t\ at fixed energies and Q 2 . 
Experimental data from DESY-HERA [T| (see text). 



FIG. 2: (Color online) The total cross section for DVCS on 
a proton as a function Q 2 at fixed energy W = 82 GeV. 
Experimental data from DESY-HERA Jrg] (see text). 



be accounted for, in the limit that x' <C x <C 1, by mul- 
tiplying the elastic differential cross section by a factor 



R g , given by [§] 



R g (X) 



2 2A+3 r(A + 5/2) 
T(A + 4) : 



with 



_ dla[A(x 

A 



\\t\)] 



dhx{l/x) 



For phenomenology on dipole cross section, we take 
two saturation models which successfully describe the 
high energy DVCS data. The first one is the impact 
parameter saturation model (hereafter b-CGC). The 
idea behind it is to introduce the impact parameter de- 
pendence into the CGC model , where now the qq — p 
differential cross section is given by = 2N(x, r; b), 
where 



I _ e -^ln 4 (BrQ sat ) 



2(7 s + ^ln^) 



: rQ sa , t < 

'■ r Qsat > 

(8) 



where Y = ln(l/a;) and Q sat (a;, 6) is the impact parame- 
ter saturation scale defined as flQl: 



tM) = ( 



x) 



exp 



-4A 

2Scgc/ 



(9) 



necessity of considerations about the impact parameter 
details of the process. The authors expressed the exclu- 
sive production cross-sections in the high-energy limit in 
terms of the Fourier-transformed dipole scattering am- 
plitude off the proton A q q{x, r, A). An important result 
about the growth of the dipole amplitude towards the sat- 
uration regime is the geometric scaling regime [l4l [l5| . 
It first appeared in the context of the proton structure 
function, which involves the dipole scattering amplitude 
at zero momentum transfer. At small values of x, in- 
stead of being a function of both the variables r and x, 
the dipole scattering amplitude is actually a function of 
the single variable r 2 Q 2 &t (x) up to inverse dipole sizes sig- 
nificantly larger than the saturation scale Q S at{x). More 
precisely, one can write 



In the b-CGC model the evolution effects are included 
via an approximate solution to the Balitsky-Kovchegov 
equation [f 21 ]. The parameters of the model can be ob- 
tained in Table 4 of Ref. [T(|. In particular, one has the 
slope Bqgc — 5.5 GeV~ 2 . The constants A and B are 
obtained from continuity conditions at rQ sat = 2. 

The second saturation model that we will consider is 
that one presented in Ref. iI3| (hereafter MPS model). 
It has the great advantage of giving directly the t de- 
pendence of elastic differential cross section without the 



A qQ (x, r, A = 0) = 2nR 2 p X {r 2 Q 2 sat {x)) , 



(10) 



2 implying the geometric scaling of the total cross-section 
at small x, i.e. a^ ot p ^ X (x, Q 2 ) = ^J^ X {t p = Q 2 /Q 2 at ). 
This geometric scaling property can be extended to the 
case of non zero momentum transfer [l6| , provided r A <C 
1. Therefore, it has been obtained that equation (fTOjl can 
be generalized to 

A qq (x,r,A)=2irR 2 p F(A)N(r 2 Q 2 sat (x,A)) , (II) 

with the asymptotic behaviors Q 2 &t (x, A) ~ 
max(Qg, A 2 ) exp[— A ln(a;)] and an unknown form 
factor -F(A) of non-perturbative origin. Specifically, the 
t dependence of the saturation scale is parametrised as 



(x,|*|) = Q 2 (l + c|t|) 



(12) 



in order to interpolate smoothly between the small and 
intermediate transfer regions. The form factor F(A) — 
exp(— B\t\) catches the transfer dependence of the proton 
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vertex, which is factorised from the projectile vertices and 
does not spoil the geometric scaling properties. Finally, 
the scaling function N is obtained from the forward sat- 
uration model [hO ] , whose functional form is given by Eq. 
([8]). The parameters for N are taken from the original 
CGC model [TTJ] . The final expression of MPS saturation 
model is given by 

Atf{x,r,A) =2wI^e- B ^N(rQ Ba . t (x,\t\),x), (13) 

which is an extension of the forward model [ll| includ- 
ing the QCD predictions for non zero momentum trans- 
fer, which reproduces the initial model for \t\ — and 
ensures that the saturation scale has the correct asymp- 
totic behaviors. The remaining parameters are c = 3.776 
GeV" 2 and B = 3.740 GeV~ 2 (obtained from a fit with 
boosted-Gaussian meson wavefunction to describe exclu- 
sive vector meson production). In order to compute the 
DVCS cross section, the amplitude above can be directly 
replaced in first line of Eq |4} . 

In order to show the reliability of the models mentioned 
above, we compare them to the recent high energy DVCS 
measurements from DESY-HERA. In Fig. [T]the differ- 
ential cross section da/dt is shown for distinct energies 
and virtualities. The solid curves represent the results for 
MPS model whereas the dot-dashed ones represent the 
b-CGC model. Both models are in agreement with data 
[13] (only systematic errors are presented), presenting dif- 
ferent trends towards large t. This is due to the different 
implementation for the t dependence of the elementary 
elastic amplitude A q q(x, r, A). A detailed analysis about 
the analytical/numerical results of the distinct models is 
beyond the scope of the present study. Finally, in Fig. 
[2] the total cross section is presented as a function of Q 2 
at fixed energy W — 82 GeV, including the experimental 
results 17]. The agreement is once again is very good. In 
this case, the numerical results for both models are ob- 
tained by integrating the differential cross section da/dt 
for \t\ < 1, which is the experimental cut. 



III. PHENOMENOLOGICAL RESULTS FOR 
ELECTROWEAK DVCS 

In this section we will extend the study for the case 
of electroweak DVCS, focusing first on the neutral weak 
current. This can occur for instance in the neutrino scat- 
tering off an unpolarized proton target through the ex- 
change of Z° instead a photon as in the electromagnetic 
case. In addition, the virtual exchange can be a with 
the production of an energetic photon, a n^, with either 
a recoiling nucleon or nucleon resonance. Using the color 
dipole picture, the elastic cross section for the high en- 
ergy neutral weak DVCS scattering on a proton target 
can be written as 

A z * p ^ p (x,Q,A) = f d 2 r [ dz (%V z o) T A qq (x,r,A) 
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FIG. 3: (Color online) The differential cross section da/dt for 
the weak neutral DVCS process (see text). 



where now {^*^z°)t denotes the overlap of the incom- 
ing Z° and outgoing real photon final state wave func- 
tions. Explicit expressions for the boson wavefunctions 
can be found in Ref. [HI, [n|. Corrections for skew- 
ness is implicitly understood, as done in previous sec- 
tion. The present formalism has been successfully used 
in calculations of neutral/charged current neutrino cross 
sectio ns J20I 2ll I22I ] and neutrino-nucleon structure func- 
1 123 1 . For the charged current case, the overlap 



= i Jd 2 r J^dz J d 2 b (*;* z o) T e 



ib-A & a qq 
d 2 b ' 



(14) 



tions If 

of the wave functions (^^z ) is replaced by w ±), 
with a possible correction for the corresponding produc- 
tion of nucleon resonance. 

We are now ready to compute the numerical results. 
For our purpose here, the real part of amplitude will be 
disregarded. For not so large Q 2 , the effective power 
X(Q 2 ) w 1.26 and the real part correction 77 ~ tan(7rA/2) 
contributes just a few percents. In Fig. [3] the differ- 
ential cross section da/dt is shown as a function of \t\ 
using expression on Eq. (|14p . In order to investigate 
the sensitivity on energy and virtuality, tree representa- 
tive cases are considered. We start with a low energy, 
W = 20 GeV and small Q 2 = 2 GeV 2 , which is suit- 
able for the neutrino beam facilities. The intermediate 
case is compatible with the DESY HERA kinematical 
range for DVCS. Finally, we extrapolate the results for 
very higher energy W = 500 GeV and large Q 2 = 50 
GeV 2 . At small \t\ < 0.4 the differential cross section 
is consistent with an exponential behavior and a Q 2 de- 
pendent slope. However, the general \t\ dependence is 
determined by the dipole cross section, da qq /dt. We 
checked that the cross section can be parameterized as 
da/dt = a (1 + c \t\) exp(— B\t\), with the constants a, c 
and B depending on Q 2 . 

In the following we compute the total cross section, 
performing the integration over |i| < 1 GeV 2 . In Fig. 
3] the energy dependence is presented for Wz P > 10 
GeV at fixed virtualities Q 2 = 1-20 GeV 2 . Roughly 
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FIG. 4: (Color online) The total cross section for weak neutral 
DVCS as a function of energy at fixed virtualities (see text). 



FIG. 5: (Color online) The total cross section for weak neutral 
DVCS as a function of Q 2 at fixed energy (see text). 



components, assumes the form [U, @, @], 



speaking, this interval probes values of x going from 
x 1CP 1 to x » 1CP 6 . The energy power of the cross 
section, a to t oc W a , ranges on a = 0.38 — 0.65, which 
increases with Q 2 . For instance, in high energies one has 
a tot {Z*p -» 7 p) = 1.55/ib(W^/GeV) - 38 at Q 2 = 1 
GeV 2 and a tot (Z*p -> jp) = 9 nb (W Zp /GeV)°- 65 at 
Q 2 = 20 GeV 2 . Therefore, the cross sections are large 
and a factor 10 2 higher compared to the DVCS results 
measured at HERA. As a short comment, the effective 
power a decreases towards low virtualities as a conse- 
quence of the non-linear effects in the dipole cross section 
since the saturation scale, Q sat (x) = Qo exp[— 0.5A ln(x)] , 
is of order 1 GeV or even larger in the kinematical range 
considered here. 

In Fig. [5] the total cross section is shown as a function 
of Q 2 at fixed energies, 40 GeV < W Zp < 500 GeV. 
The virtuality dependence is similar to that appearing 
on DVCS results at HERA. Once again, the absolute 
values are large at high energies reaching to atot — 10 
fib for Q 2 = 1 GeV 2 . In the interval 1 < Q 2 < 10 
GeV 2 the cross section can be parameterized as a to t — 
a/(Q 2 ) s , with the constants a and S s» 1 depending on 
W. When going out of this interval, the Q 2 dependence 
is steeper and more involved. For example, one obtains 
a to t{Z*p -> tp) = ^fjb/iQl/Q 2 ) 11 for W Zp = 40 GeV 
and a to t(Z*p -> 1P ) = 11 fib (Q 2 /Q 2 ) 97 at W Zp = 200 
GeV, where Q 2 , = 1 GeV 2 . 

The present calculation can be serve as an input for 
electroweak DVCS in neutrino interactions. In general, 
a neutrino induced DVCS process in a nucleon is given 
by the reaction h (k) + N\ (pi) -> h{k') + N 2 (p 2 ) +7(92), 
where a neutrino £1 scatters from a nucleon N\ to a final 
state h, nucleon N 2 and real photon 7. The differen- 
tial cross section in the target rest frame, in which the 
weak virtual boson four-momentum q has no transverse 



d 4 a 



l + x(M/E u ) 



dxdQ 2 d\t\dct> 64 s (2tt) 4 MEuX [y + 2x(M/E v )Y 



in 



where T cx A(x, r, A) represents the invariant matrix el- 
ement. The invariants in equation above are given by 
s = (k + Pl ) 2 , Q 2 = -q 2 and y=(pi- q)/(pi ■ k). The 
quantity E v denotes the energy of the incoming neutrino 
beam and (A the angle between the lepton and nucleon 
scattering planes. In electromagnetic DVCS, T receives 
a large contribution from the Bethc-Hcitlcr (BH) process. 
Considering a neutrino scattering off an unpolarized pro- 
ton target through the exchange of Z n , one only measures 
the Compton contribution since the photon can not be 
emmited by the neutrino. In the case of weak charged 
current case, where neutrino-nucleon scattering occurs 
through the exchange of W + with a proton in the final 
state, one has both DVCS and BH contributions since 
the outgoing muon can emit the real photon. 

Finally, we will present a phenomenological investiga- 
tion on the weak DVCS process for neutrino scattering 
off nucleus. In order to do so, we will rely on the geo- 
metric scaling property [TJ, [l5[ of the saturation models 
within the color dipole approach. The geometric scal- 
ing means that the total j*p cross section at large ener- 
gies is not a function of the two independent variables 
x and Q, but is rather a function of the single vari- 
able t p = Q 2 /Q 2 at (x) as shown in Ref. Q. That is, 
<Jy* p {x,Q 2 ) = <t 7 * p (t p ). In Refs. [l5| it was shown that 
the geometric scaling observed in experimental data can 
be understood theoretically in the context of non-linear 
QCD evolution with fixed and running coupling. Re- 
cently, the high energy l^p, pA and AA collisions have 
been related through geometric scaling 24j. Within the 
color dipole picture and making use of a rescaling of the 
impact parameter of the "/*h cross section in terms of 
hadronic target radius Rh, the nuclear dependence of 
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the j* A cross section is absorbed in the ^-dependence 
of the saturation scale via geometric scaling. The rela- 
tion reads as a] ot (x,Q 2 ) = ka clot (<3sat, P -> Q sa t,A), 
where ka = {Ra/ Rp) 2 ■ The nuclear saturation scale was 
assumed to rise with the quotient of the transverse par- 
ton densities to the power A « 1 and Ra is the nuclear 
radius, Q 2 atA — (A/ka) a Q 2 at , p - This assumption suc- 
cessfully describes small- x data for ep and eA scattering 
using A = 1.26 and a same scaling curve for the proton 
and nucleus (24|. 

Following the simple arguments of scaling, we replace 
Rp -> R A in Eq. (T3|) and QL, P -> (AR 2 p /R 2 A ) A Q 2 sat , p 
in Eq. (TT2]). In case of A = 1 the previous replace- 
ment becomes the usual assumption for the nuclear sat- 
uration scale, Q 2 atA — Qsat, P - The corresponding 
phenomenological results are presented in Fig. [5] As 
a short comment, the present estimation corresponds to 
the coherent electroweak DVCS process, Z°A — > jAX, 
where the nucleus remains intact after interaction [26| . 
Future studies should include the incoherent contribu- 
tion, which is sizeable for high energies and heavy nuclei 
[26| . The total cross section as a function of energy is 
shown for fixed virtuality Q 2 = 5 GeV 2 and representa- 
tive nuclei. The effective power on energy is dependent 
on the mass number. For light nuclei like carbon the 
power follows similar trend as for the nucleon case but 
there is a clear suppression at high energies. For large nu- 
cleus, i.e. Ca, Fe or Pb, the effective power systematically 
diminishes due to nuclear shadowing. In the saturation 
model, the enhancement of the nuclear saturation scale 
due to the presence of nuclear targets leads to an addi- 
tional flattening of the energy growth. The nuclear sup- 
pression as a function of energy is verified by contrasting 
the nuclear results to the proton case (upper dot-dashed 
curve). In order to quantify the size of nuclear suppres- 
sion, let us give some examples. For carbon it is of order 
15% at high energies, whereas for lead it reaches 50 %. 
The present study for nuclear targets should be ti mely 
for future neutrino- nucleus experiments as Minerva |25j ]. 

IV. SUMMARY 

We estimate the differential du/dt and total cross sec- 
tion for the high energy electroweak deeply virtual Comp- 
ton scattering on a proton target and extended for nu- 



clear targets. We take into account the QCD color dipole 
formalism, which is quite successful in describing small- a: 
data. The absolute values for the cross sections are large, 
higher than electromagnetic DVCS measured at DESY- 
HERA, which makes its measurement feasible. The main 
theoretical uncertainty is the specific model for the dipole 
cross section. We have considered two saturation model, 
which are able to correctly describe electroweak DVCS 
data. Skewedness correction have been properly included 
in calculations, giving a sizable enhancement. In the 
present study, one considers neutrino scattering off an 




W [GeV] 

FIG. 6: (Color online) The nuclear cross section for weak 
neutral DVCS as a function energy at fixed Q 2 = 5 GeV 2 . 
Results for representative nuclei are presented, (see text). 



that is neutral current case. The extension for charged 
weak current is possible within the dipole approach. 
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